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Analyzing the m_|_-spectrum and two-particle correlations of negative hadrons from 



158 AGeV/c Pb+Pb collisions at slightly forward rapidities we find a (thermal) freeze-out 
temperature of about 100 MeV and transverse flow with v± fh 0.55c. The Mj_-dependence 
■ of the correlation radii prefers a box-like transverse density profile over a Gaussian. From 
an analysis of the pion phase-space density we find fj, n ~ 60 MeV at thermal freeze-out. 

1. MOTIVATION 

In relativistic heavy-ion collisions the momenta of the escaping hadrons are fixed at 
the point of "thermal freeze-out". They carry information about the final state of the 
fireball which one can use to reconstruct this state and thus obtain boundary conditions 
for back-extrapolation into earlier stages of its evolution. In this way one can examine 
^ , the question whether or not a quark-gluon plasma was created. 

Assuming thermalization and collective transverse flow at freeze-out, these two features 
can be extracted separately by combining single-particle m^-spectra and two-particle 
Bose-Einstein correlations PJfl. Particles of non-zero transverse momentum are emitted 
from a fireball region which moves outwards transversely; in the local comoving frame their 
transverse momentum is lower, resulting in an enhanced Boltzmann factor. Transverse 
flow thus flattens the m^-spectrum. A given spectral slope can, however, be reproduced by 
different combinations of flow and temperature, because larger flow can be compensated 
for by lower temperature and vice versa. On the other hand, the sizes of the regions 
emitting particles of given momentum (homogeneity regions, measured by correlation 
radii) are controlled by the expansion velocity gradients. The natural velocity measure 
is given by the chaotic thermal motion; increasing T and the transverse flow (velocity 
gradients) together at the proper rate will thus not change the resulting correlation radii. 
In summary, T and v± are anticorrelated by a given spectral slope, but correlated by the 
measured homogeneity sizes. Exploiting both types of correlations their values can be 
determined unambiguously. This method was used for analyzing the data from central 
Pb+Pb collisions at 158 AGeV/c taken by the NA49 collaboration gj. 
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2. THE MODEL 

It is impossible to extract from Bose-Einstein correlations the global fireball features in 
a model-independent way. We therefore parameterize the fireball by an emission function, 
i.e. by a model for the (Wigner) phase-space density of the particles at freeze-out |l|J^,|]: 

4 m± cosh(?/ - 77) / K-u(x)-/i\ ( {v - Vo) 2 \ 
b(x,K)a x = — — exp — exp 



(2tt) 3 " \ T J " V 2 ( A v) 
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xG(r) j== exp ("T^J (1) 

The main assumptions here are the following: (i) Approximate thermal equilibrium at 
temperature T. (ii) Boost-invariant longitudinal flow combined with a tunable transverse 
expansion profile with a transverse flow rapidity given by i] t (r) = r]f(r/r ims ); this expan- 
sion is encoded in the velocity field u(x). (iii) For the transverse density profile G(r) two 
different forms were studied: 

( r M 

box-like: G(r) = 9(Rb — r) , Gaussian: G(r) = exp —pr . (2) 

V 2R gJ 

If the parameters T, rjf, and Rb or Rq are independent of space-time rapidity rj, the 
model cannot reproduce the rapidity dependence of the data. We therefore concentrate 
here on a single rapidity window (3.9 < Y < 4.4). The time parameters r , At were fixed 
together with the other parameters from the correlation radii, whereas the longitudinal 
width Arj = 1.3 of the density profile was fixed by adjusting the width of dN/dy. 

The correlation radii were calculated numerically from the "model-independent expres- 
sions" p|. When calculating the single-particle spectrum resonances were taken into 
account as in [0; since we fitted the h~ spectrum we also added negative kaons and an- 
tiprotons. For the resonances baryonic and strangeness chemical potentials, /zb and //g, 
were taken into account as required by baryon number and strangeness conservation |J. 

3. THE FIT 

Fig. [I] shows the x 2 contours of our fit to the data. The valley extending from the 
upper left to the lower right corner shows the ambiguity of the extraction of temperature 
and transverse flow from the spectrum. Moreover, the left and right figure show also a 
model ambiguity: at fixed temperature different density profiles lead to slightly different 
transverse expansion velocities. The results of the fit to the correlation radii (width 
parameters in a Gaussian parametrization of the correlation function) indicate that a 
box-like density fits the data better than a Gaussian one, although the latter cannot 
be excluded. Strictly speaking, only models with too weak transverse flow seem to be 
completely wrong. With the box-model we find T w 80 — 110 MeV and v± m 0.47c — 0.62c. 
The radius of the transverse density distribution Rb = 12.1 fm, corresponding to r rms = 
8.6 fm, is about twice as large as the initial transverse overlap region of the two Pb nuclei, 
consistent with the strong transverse flow. 

Both models fit the m^-spectrum well but, as shown in Fig. ^|, they differ in the quality 
of the fit to the correlation radii. The M^-dependence of transverse correlation radii R Q 
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Figure 1. \ 2 contours resulting from the fits to the single-particle m_|_-spectrum of h~ 
|| (dashed) and the two-particle correlations (solid). The contours correspond to the 
indicated confidence levels. The two plots are obtained with different transverse density 
profiles. Crosses show the positions of parameter sets used in comparative calculations. 

and R s is better captured by the box-model. This can be traced back to the fact that the 
box-model has a rigid radial density cutoff; for increasing K± the homogeneity region is 
thus squeezed towards that edge, making it narrower in the direction of K±. As a result 
R Q decreases more strongly with K± than for the Gaussian model where the homogeneity 
region expands into the tail of the density distribution (see Fig. ||b). We stress that it is 
the generic feature shown in the upper half of Fig. 0b (and not the box-like transverse 
density distribution itself) which appears to be required by the data. 



4. MULTIPLICITY 



From the fitted models we can calculate the total pion multiplicity, assuming chemical 
equilibrium at thermal freeze-out with /iB ~ 360 MeV, corresponding to an entropy per 
baryon of about S/A=38 @,(|, and /is fixed by strangeness neutrality. This implies fi n =0 
and results in a predicted multiplicity which is more than a factor 3 below the measured 
value. By investigating the pion phase-space density || we inferred a chemical potential 
for the direct pions of fj, n ~ 60 MeV at T « 100 MeV. This implies that at thermal freeze- 
out the pions are out of chemical equilibrium. This is expected if chemical freeze-out 
happens in equilibrium near the confinement phase transition at T CT « 170 MeV [pLIlO . 



5. CONCLUSIONS 



A combined analysis of single-particle spectra and two-particle correlations avoids the 
ambiguities and strongly reduces the model-dependencies resulting from an analysis of the 
m_i_-spectrum of only a single particle species flD] . It leads to the unavoidable conclusion 
that thermal freeze-out of the fireball happens at a very low temperature (T t herm ~ 80 — 
110 MeV) and is accompanied by very strong transverse flow (v± ~ 0.47 c — 0.62 c). The 
M^-dependence of the homogeneity lengths can be somewhat better reproduced by a 
box-like than by a Gaussian transverse density profile at freeze-out. 
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Figure 2. (a) Comparison of the measured correlation radii with the two model fits. The 
temperature (T=100 MeV) and transverse flow (v± = 0.5 c for box, v± = 0.46 c for Gauss) 
of the fits are marked by crosses in Fig. [I], (b) Transverse cuts of the homogeneity regions 
in the two models for pions with K± = 400 MeV. 



At this low temperature the fireball is no longer in chemical equilibrium. The pion 
abundance appears to be fixed by chemical freeze-out already at T chem « T cr « 170 MeV 
[@,[l(J. We find that at T therm =100 MeV this leads to a pion chemical potential of fi n ~ 60 
MeV §JTJ. 
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